Localized transverse bursts in inclined layer convection 
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We investigate a novel bursting state in inclined layer thermal convection in which convection rolls exhibit 
intermittent, localized, transverse bursts. With increasing temperature difference, the bursts increase in duration 
and number while exhibiting a characteristic wavenumber, magnitude, and size. We propose a mechanism which 
describes the duration of the observed bursting intervals and compare our results to bursting processes in other 
systems. 

PACS numbers: 47.54.+r, 47.20.-k, 



Bursting phenomena are a common feature of nonequilib- 
rium systems. Various bursting mechanisms have been de- 
scribed by Knobloch and Moehlis [[!]] and characterized in 
terms of their recurrence properties and dynamic range, distin- 
guishing between temporally localized bursts which are spa- 
tially global and those which are spatially local. Some ex- 
amples of the former are the breakdown of spiral vortices in 
Taylor-Couette flow P, pfl, fluctuations in heat transport in bi- 
nary fluid convection [^T||], and shear flow turbulence |^, [7]]. 
While a dynamical systems approach has been fruitful in the 
global case, spatially localized bursts remain less understood. 

In the Letter, we characterize a recently discovered [|[| 
bursting phenomenon in inclined layer convection (ILC) re- 
sulting from the presence of both buoyancy and shear instabil- 
ities. Two features are intriguing: spatially localized bursts of 
a characteristic size and the triggering of local disorder within 
the burst without completely destroying the underlying roll 
structure (see Fig. [j]). Temporally, we observe multiple cy- 
cles of transverse modulation, turbulence, and decay within 
the bursts. Existing theory does not address either the spatial 
or temporal dynamics. 

Inclined layer convection is a variant of Rayleigh-Benard 
(thermal) convection [|T(J, in which a fluid layer is heated 




FIG. 1: (a-d) Sequence of background-divided shadowgraph images 
and associated contrast-enhanced power spectra for Cell 1 at 9 = 78° 
and e = 0.08. (a) quiescent, (b) modulations during early bursting, 
(c) developed bursting, and (d) end of burstlet. Times correspond to 
points marked in Fig. n. Uphill direction is at left, (e) Disordered 
bursting at 6 = 78° and e = 0.09. Also see movies at [§]. 




FIG. 2: Schematic diagram of inclined layer convection profile. 

from one side and and cooled from the other. Tilting this 
layer by an angle 9 (see Fig. |]) results in a base state which 
is a superposition of a linear temperature gradient and a shear 
flow up along the hot plate and down along the cold. When 
heated beyond a critical temperature difference AT C , the fluid 
convects due to the buoyancy of the hot fluid. As 9 is in- 
creased, the buoyancy provided by the perpendicular compo- 
nent of gravity g± = g cos 9 becomes weaker and the shear 
flow provided by becomes stronger. Above a codimension- 
two point at 9 C , the primary instability is due to this shear flow 
instead of buoyancy [[TTj [l2[ OJ. Interesting bursting behav- 
ior has been observed in the vicinity of this codimension-two 
point [§. 

We perform experiments in high pressure CO2 in an ap- 
paratus similar to that described in modified to allow 
for inclination. The gas was at a pressure of (48.26 ± 0.01) 
bar regulated to ±0.005 bar with a mean temperature of 
(25.00 ± 0.05)°C regulated to ±0.3mK. Our three convec- 
tion cells were of height d = 778 ± 2/im and length 97.3<i. 
The widths in the x direction are Li = 10. Ad for Cell 1, 
L 2 = 20.9d for Cell 2, and L 3 = 3l.0d for Cell 3. These pa- 
rameters give a vertical viscous diffusion time of t v — d 2 jv = 
(4.50 ± 0.02) sec, a Prandtl number a = 1.301 ± 0.001, 
and weakly non-Boussinesq conditions (Q = 0.2 to 0.8, as 
described in Jiot l for horizontal fluid layers). The planform 
of the convection pattern was observed via the shadowgraph 
technique Jl4} | using a digital camera. The two control param- 
eters are the angle 9 and the nondimensionalized temperature 
difference e = 1- Images were collected at 27 frames/r„ 

with runs of duration at least 1000 t v at various values of e, 9. 

Fig. U shows the phase diagram for ILC in the vicinity of the 
codimension-two point, with onset measured experimentally 
(data points) and compared to predictions from linear stability 
analysis JO}], The localized transverse bursts appear intermit- 
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FIG. 3: Phase diagram near the codimension-two point of ILC, for 
Prandtl number a — 1.301. Stability curves are from ; solid lines 
are primary instabilities and dashed lines are extensions. Shaded area 
corresponds to data collection on bursts. 
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FIG. 4: Wavenumber of transverse bursts p (determined from power 
spectra) as a function of 8, e, and L. In each graph the data plotted 
are averaged over the other two variables. Dashed lines are p c /2 for 



tently within the underlying longitudinal rolls as a secondary 
instability for 8 < 9 C . These bursts are triggered within a 
secondary instability in which there is already local growth of 
regions of high-amplitude longitudinal convection, as seen in 
the sequence of images in Figju a movie of the corresponding 
images are available online [gj. Within these regions, trans- 
verse modulations repeatedly grow and decay, becoming dis- 
ordered or turbulent in the process. Eventually, this cycle of 
bursting ends and the system returns to quiescent, weak lon- 
gitudinal rolls. The phenomenon becomes more pronounced 
at large e, so that eventually the whole cell is bursting. Each 
localized patch does not spread and no global bursting was 
observed. At lower inclination (further from the codimension- 
two point), the bursting is indistinguishable from the crawling 
rolls described in [^J . 

The measured wavenumber p of the transverse modulations 
is shown in Fig. ^| and observed to be approximately con- 
stant. At — 80°, the transverse rolls were determined to 
have a wavenumber of p c = (2.82 ± 0.04)rf at onset; the 
burst wavenumber p is nearly subharmonic with respect to this 
value. The appearance of this bursting phenomena close to 
the stability curve for transverse rolls suggests that a related 
shear instability is playing a role. Nonetheless, no mode reso- 
nant with p c is observed as part of the bursting spectrum (see 
Fig.[|. 
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FIG. 5: Burst width in x (gray) and y (black) directions: (a) vs. L at 
onset of bursting for 6 = 78°, (b) vs. 6 at onset of bursting in Cell 2, 
and (c) vs. e in Cell 2 at 6 = 78° . 



Using Fourier decomposition, we determined the area of 
each cell occupied by transverse bursting and found the width 
of the region for those instances in which only a single burst 
was present. As shown in Fig. |5[ the bursts generally have a 
characteristic of size of 15d in the transverse direction and 20d 
in the longitudinal. Bursts are smaller in Cell 1 (see Fig. ||a), 
where they fill the cell in the x (transverse) direction, and 
larger for 9 — 76° (Fig. |^b) and increased e (Fig. ||c) Cell 
2 and Cell 3 typically contain multiple bursts; therefore, we 
will utilize data from Cell 1 to focus on the temporal behavior 
of single bursts. For Cells 1, 2, and 3 the onset of bursting 
was observed to occur at e = 0.075, e = 0.02, and e = 0.04, 
respectively [pjl]. 

We describe the bursts based on the modes they excite in 
Fourier space, using the power spectra of background-divided 
images as a function of time. Sample spectra and their corre- 
sponding shadowgraph images are shown in Fig. [j], with three 
prominent modes: a pure x mode q, a pure y mode p, and 
a mixed mode pq. The longitudinal rolls are composed of the 
pure q mode, while the transverse modulations are constructed 
from the p and pq modes. Higher-order modes serve to cre- 
ate the characteristic shape of the bursts and will be ignored 
here. A further simplification should be noted: the shadow- 
graph technique integrates through the fluid layer (z), while 
these bursting phenomena exhibit three dimensional behavior, 
particularly when the modulations become disordered. 

We define S q , S pq , and S p as the power in each of these 
three peaks, taken as the total intensity in a fixed region of 
Fourier space. Fig. || shows time traces of the power in each 
of the three modes. Prior to the beginning of the burst S q typi- 
cally increases in relative amplitude, as shown in Fig. [7|a. The 
beginning of a burst is characterized by the rapid growth of 
the transverse modulations (via S pq ). The S p mode is weaker 
and grows after a delay of about 0.5t v (see Fig. |7|b), which 
suggests that it is due to a nonlinear interaction or resonance 
condition between the q and pq modes. The simplest such 
nonlinearity would be quadratic, and by multiplying S q (t) and 
S pq (t), we obtain a good reproduction of the curve for S p (t), 
shown as the gray dashed line in Fig. ^. 

The time series in Fig. || show two temporal features: (1) 
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FIG. 7: (a) Cross-correlation of (S pq , S q ), dashed and (S pq , S p ), 
solid, (b) Ensemble average (S(t — t s t ar t)) over all bursts where 
t&tart is the rising edge of S pq for each burst. Both are plotted for 
Cell 1 at e = 0.08 and 9 — 78°. Power is in same arbitrary units as 
Fig. ft 



intermittent, alternating periods of quiescence and bursting 
and (2) burstlets characterized by peaks and troughs within 
the bursting intervals. As S pq (t) rises at the beginning of 
a burstlet, the pattern exhibits the increasingly well-defined 
modulations shown in Fig. |l]b. When S pq peaks, these modu- 
lations become spatially disordered on short time scales, mov- 
ing rapidly within the underlying rolls. While the general roll 
pattern is retained in any single image (see Fig. [l|c), the roll 
segments move turbulently within the localized bursting re- 
gion [Pp. At higher e and 9, both of which increase the shear 
flow, this disorder is increased (see Fig. |l]d). From the disor- 
dered state, the modulations either decay — resulting in the 
end of the bursting interval — or grow again, creating another 
burstlet within the same interval. A succession of such events 
is suggestive of the presence of a limit cycle. Stochastic limit 
cycles have been previously described in [|IJ |l7|, relying on a 
random forcing effect such as pressure fluctuations to produce 
random events with a single, well-defined mean rate. 

To determine the quiescent and bursting intervals, we set a 
threshold for S pq (t) above which the system was considered 
to be bursting. Because the onset and decay of the burst is 
sharp, the results were not sensitive to the choice of thresh- 
old over a reasonable range of values. Using this information, 
we examined the duration of quiescent intervals Tq, bursting 
intervals Tb- Burstlets are separated by intervals t\, within 
the Tb, with burstlet peaks identified by local maximum. Ex- 
amples of the determined bursting intervals are shown by the 



gray bars at the bottom of Fig. |[ Sample probability distri- 
bution functions (PDFs) and mean values as a function of e 
are plotted in Fig. ^| The quiescent intervals Tq were longest 
close to the onset of bursting, and decreased until they were 
no longer detectable at higher e. Conversely, the bursting in- 
tervals Tb grow with e. The combined effect is that at high 
e the whole cell is in a perpetual state of bursting since each 
localized burst lives longer and new ones begin sooner. All of 
these trends hold at other values of e and 9 as well. 

Such behavior distinguishes this bursting from the behav- 
ior of bursts in Taylor-Couette flow |Q, |[j . There, the bursting 
is attributed to a secondary instability whose growth above a 
threshold triggers turbulence throughout the fluid. Once the 
turbulence has begun it destroys the underlying rolls which 
were providing energy and thus it dies away after a well- 
defined period of time. It appears that the initial transverse 
modulations in ILC are such a secondary instability, although 
they appear to arise due to a growing primary mode which 
triggers their onset, as shown in Fig. 0. However, instead of 
fully bursting, the burstlets cause only local disorder and fail 
to globally trigger turbulence. Because the modulation and 
rolls are not completely destroyed (see Fig. ft), they are read- 
ily able to grow back up again after they decay. This creates 
bursts comprised of multiple burstlets. As a result, Tb is not 
constant as was observed for the Taylor-Couette bursts but in- 
stead increases with e, as shown in Fig. ftf. The mechanism 
for the growth of the Taylor-Couette bursts [||, |[| provides for 
(Tq) ~ 1/e behavior based on a constant growth rate of the 
secondary instability. It is possible that a similar mechanism 
is at work since a similar trend of longer quiescent periods at 
lower e is observed (see Fig. 

To create cycles of burstlets within each bursting interval, a 
second bursting process must be at work as well. At each e and 
9 the PDF of T(, (see Fig. [|b) exhibits a negative exponential 
distribution 
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where (rt) is the mean waiting period of a Poisson process 
generating the bursts. Because the lowest Tb bin is underre- 
ported due to the finite sampling of the time series, (Tb) e ff 
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FIG. 8: (a,b,c) Sample PDFs of Tq, Tb, and N respectively at e = 
0.08 and = 78°. (b) Dashed line is exponential distribution fit to 
data excluding the first bin. (c) Dashed line is Eqn. ^ plotted from 
the measured mean, (d) Mean quiescent period as a function of e at 
9 = 78° . (e) Measured and effective (from fit) mean burstlet period 
and (f) mean number of burstlets and mean bursting interval as a 
function of e at 9 = 78° . 



is found instead by fitting the negative exponential to all data 
except for the first bin. 

This description can also explain why Tb isn't constant. 
In a Poisson process, the events are memoryless: the wait 
time for the next event is independent of how long the sys- 
tem has already been waiting. If a new burstlet were not 
generated before the decay to quiescence, the modulations 
would die away and the burst interval would be over. As- 
suming a constant, unknown decay time tj, the probability 
distribution for N burstlets for which each T{, is no more than 
Td is calculated using the cumulative distribution of Eqn. [j] 
V(n < T d ) = 1 - e- T »/ Td . 



P(N) = 



(1 - e -1 "*/^)^ 

e r d /(r b ) _ I 
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This distribution contains a single parameter (N) = e Td l ^ Th \ 
which can be measured directly from the experimental V(N). 
Eqn. ^ then reduces to 
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terval. This formulation allows for comparison with the ob- 
served data with no fit parameters. Fig. §b shows good agree- 
ment with this model, as do plots at other parameter values. 
Therefore, our results are consistent with the burstlets being 
Poisson-distributed events. 

Transverse bursts in inclined layer convection show many 
consistent and unexplained properties across a range of e and 
6: spatial localization, transverse modulations of uniform 
wavenumber, persistence of the roll structure while bursting, 
and repeated cycles of growth and decay once triggered. As- 
pects of both spatiotemporal chaos and turbulence appear to 
be relevant and these results represent a first step in under- 
standing the structure and dynamics of this novel bursting 
phenomenon. In particular, more work is needed to under- 
stand both the local growth of longitudinal rolls which triggers 
the bursting and the implications of a resonance condition be- 
tween the p, q, and pq modes. 
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where (N) is the measured mean number of burstlets per in- 
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